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Abstract —In this paper, we propose low complexity oppor¬ 
tunistic methods for Interference alignment in -transmitter 
MIMO Interference channels hy exploiting multiuser diversity. 
We do not assume avallahility of channel state information (CSI) 
at the transmitters. Receivers are required to feed hack analog 
values Indicating the extent to which the received interference 
subspaces are aligned. The proposed opportunistic interference 
alignment (OIA) achieves sum-rate comparable to conventional 
OIA schemes but with a significantly reduced computational 
complexity. 

Index Terms —Interference alignment, user selection, user pair¬ 
ing, sum-rate, multiple-input multiple-output (MIMO), interfer¬ 
ence channel 

I. Introduction 

Interference alignment (lA) is a promising interference 
management technique for future wireless networks which 
are interference limited, such as, MIMO interference channels 
(IC), MIMO interfering broadcast channels (lEBC), etc. It was 
demonstrated in HI that lA can achieve a sum degree-of- 
freedom (DoE) of -L in a AT-user SISO interference channel. 
lA utilizes multiple signaling dimensions (due to multiple an¬ 
tennas or time/frequency extensions) to suppress the received 
interference into a reduced dimensional subspace of the receive 
space. 

Conventional methods for interference alignment Q-il de¬ 
pend on one or more of global channel state information, chan¬ 
nel state information at transmitters, reciprocity of downlink- 
uplink channels, transmitter cooperation or iterative methods. 
If these assumptions are relaxed, it is not possible to achieve 
interference alignment by employing transmit precoding in 
order to align the interferences received at the receivers and 
thus we rely on opportunistic methods to select users for which 
the interferences are naturally aligned. 

Low overhead feed back based OIA has been proposed 
in Q-Jll. A 2 X 2 MIMO IC with 3-transmitters has been 
considered in a, while a, a extend it to the case of 
M X 2M MIMO IC, again with 3-transmitters. In these works, 
it is assumed each transmitter has a separate user group. In 
each group, a single user is selected by the corresponding 
transmitter for opportunistic lA. In this paper, we extend OIA 
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for the general case of AT-transmitter MIMO IC. Eurther, in 
order to better exploit the available user diversity, we consider 
the problem of user pairing for achieving opportunistic lA. We 
use a geometric interpretation of the signal space to define the 
measure of alignment which quantifies the extent to which the 
interference subspaces are aligned. Each transmitter broadcasts 
a reference signal and receivers calculate their corresponding 
measure of alignment and feed it back. Depending on the 
received values for the measure of alignment, the transmitters 
can select their user independently (user selection) or they can 
be paired with the users by a central node (user pairing). 

The remainder of this paper is organized as follows. In 
Section im we describe the system model. Sectionllllldescribes 
the proposed choice for measure of alignment. Section |IV] 
discusses about OIA in the user selection framework while 
Section lY] discusses OIA in the user pairing framework. 
Performance comparison is presented in Section IVII Einally, 
Section I VII I concludes the paper. 

II. System Model 

We consider a network with AT transmitters and N receivers. 
Each transmitter is equipped with Nt antennas and each 
receiver is equipped with Nji antennas. In a AT-transmitter 
MIMO IC, each user receives (AT — 1) interfering signals and 
one desired signal, each of dimension Nt- We let Nt = M 
and Nr = 2M so that M dimensions can be designated for 
the desired data streams and the remaining M dimensions for 
interference alignment, at each user. We consider two different 
system models, namely, user selection model and user pairing 
model. 

A. User Selection 

In the user selection framework, we assume that there are 
AT cells, each with a single transmitter (base station). The 
receivers (users) are arbitrarily divided into K groups of size 
S = N/K, where N is the total number of users in the 
network. The signal £ i^Nrxi j-gceived by the nth user 
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Fig. 1. User Selection Model : Each transmitter selects and serves a single 
user in each group 

in the fcth cell is given by : 

( 1 ) 

desired signal y 

interference signals 

where , e is the channel gain matrix between 

Zth base station (BS) and user n in cell k and with each entry 
assumed to be independently and identically distributed (i.i.d.) 
circular symmetric complex Gaussian (CSCG) random vari¬ 
able with unit variance 1), w* G (je^otes an ad¬ 

ditive white Gaussian noise (AWGN) with ^ CA/'(0, Iath) 
andx; G is the signal vector transmitted by transmitter 

I, encoded by a Gaussian codebook. 

Fig. 1 depicts the user selection model. In each cell, out of a 
total of S users, only one user is selected by the corresponding 
BS for data transmission. This selection is carried out such that 
opportunistic lA can be achieved and the exact procedure will 
be discussed in Section ITVl 

B. User Pairing 

In the user pairing framework, we assume that each of the 
K BSs can connect to any user in the network. In other 
words, BSs and users are not divided into distinct cells. 
For differentiating from the user selection model, we let 
Gn,k G {Q (}jg channel gain matrix between BS 

k and user n. The received signal y„ G at receiver n 

is given by 

K 

Yn = ^ + Wn (2) 

k^l 

with each entry of G„ ^ being i.i.d. CSCG random variable 
CAf{0, 1) and w„ G is AWGN at user n with w„ ^ 

CAf{0,lN^) 


Fig. 2. User Pairing : Each transmitter is paired with one user by the center 

We define the N x K pairing matrix P with Pij as the 
entry on its ith row and jth column which is given by ; 

1, if receiver i and transmitter j form a pair 
0, otherwise 

where Pi,j — 1 Pij < 1 to ensure that each 

BS is connected to exactly one user and that each ‘connected’ 
user is connected to exactly one BS. The received signal at 
the nth user, y„ can now be decomposed as : 

K K 

Yn — 'y ^ Pn^kGji^k^k P ^ ^ (1 Pn,k)Gji ]^'K}^ (3) 

k=l k=l 

^ ^ ^ ^ ^ 
desired signal interference signals 

Fig. 2 depicts the user pairing model. The pairing is enabled 
by the presence of a central system connecting each BS 
through backhaul links. The exact procedure for finding the 
transmitter-receiver pairing configuration or equivalently P 
such that opportunistic lA can be achieved, will be discussed in 
Section |V] After the pairing configuration is found, the center 
can redirect the users’ data to their corresponding BSs. 

In both the frameworks, since the transmitters do not have 
channel state information, an equal power allocation among 
Nt data streams is assumed such that x/ ^ CJ\f{0, 
where P is the transmit power constraint, assumed common, 
at all the BSs. Note that, if transmitter I wants to convey d data 
streams where d < Nt it can employ an arbitrary precoding 
matrix V, G such that Vf V, = 1^ so that the system 

model remains statistically identical. 

III. Measure oe Alignment 

Recall that in a iT-transmitter MIMO IC, each user will 
receive {K — 1) interference signals and one desired signal, 
each of dimension Nt- In order to quantify the suitability of 



















a user for data transmission in the interference channel, we 
need to define a measure of alignment which quantifies the 
extent to which the received interference signals are aligned 
at each user. To that end, we will briefly review Grassmann 
manifold. 


A. Grassmann Manifold 

The Grassmann manifold Gn,m{'C) has been defined as 
the set of all M-dimensional subspaces of complex Euclidean 
A^-dimensional space, il. It is a widely used geometric 
concept in wireless communications and helps in the design 
and analysis of different methodologies. Let A G GN,MiC) be 
a M-dimensional subspaces. A N x M matrix A is defined 
as generator matrix for A if its columns span the subspace 
corresponding to A and forms an orthonormal bases for the 
same, i.e., A^A = 1m- Note that there can be infinitely many 
generator matrices for a given subspace A which can be ob¬ 
tained by the transformations A AU where U G C^xM 
an arbitrary unitary matrix. For subspaces A,B G Gn,m{^)^ 
the M ordered principal angles, 9i,92, ■ ■ - Om G [0,7r/2] 
between the subspaces are obtained sequentially as 

cos{9m) = max |a^b| = |a^b„| 

a.GA 

bsB 

s.t. ||a|| = ||b|| = 1 

a^a„ = b^b„ = 0 V n G Mm (4) 

where Mm = {1, 2,..., m — 1}, while {am}m=i 
{bm}m=i *^he principal vectors for A and B, respectively. 

The chordal distance between A and B is defined as 


dc{A,B) 


'i 


M 

^ sin^{9m) 

m—1 


(5) 


Alternatively, chordal distance between the two subspaces A 
and B can be represented in terms of their generator matrices 
as 


dc{A,B) = 4(A,B) = — ||AA" - BB"||f (6) 

V2 


= JM - tr(A^BB^A) (7) 


Chordal distance is known to be proportional to the degree 
of orthogonality between the subspaces. Note that, the chordal 
distance is invariant to the choice of generator matrices. 


B. Spread of Subspaces 

Let be L matrices of size N x M. Thus, 

each of these matrices would correspond to planes/subspaces 
e Gn,m{C) with {ni}ti e as their 

generator matrices. In order to define the spread of these 
subspaces, consider the following 

L 

F= argmin Vd2(F',H,) (8) 

F' s.t. F'«F'=Im 


Let T G Gn^m{C) be the plane corresponding to F. Thus, 
T can be considered as the mean of subspaces - Quite 

naturally, we can define the spread of these subspaces as 


/ = min 

F' s.t. F'»F'=Im 


E'^c(f',ho 

1=1 


The problem in ([8]l can be simplified as follows ; 


(9) 


F= argmin ^d^(F',Hi) 

F' s.t. F'"F'=Im 

L 

= argmin '^{M- Tr(F'-f^HiHf F')) 

F' s.t. F'«F'=Im 

L 

= argmax V Tr(F'^HiHf F') 

F' s.t. F'«F'=Im 

= argmax Tr f F'^ f V ) F'| 

F' s.t. F'«F'=Im V \i=i ’/ / 

= [vi(B„) V 2 (B„) ... vm(B„)] 


( 10 ) 


K _ _ 

where B = X) and v;(H) denotes the left singular 

i=i 

vector of H which corresponds to the Zth largest singular value 
ID. Thus, / in (01 is given by 

L 

f = J2dliF,Ui) 

1^1 


f = LM- Tr I^F^ j F^ 

M / L \ 

/ = LM- y] 

m^l \l^l ) 

2M / L \ 

/§ y: A„ y:H,H» 


( 11 ) 


m=M+l \l—l 


where Am(H) denotes the mth largest singular value of H 
and (i) holds because 

2M / L \ / L \ 

EA- =Tr 

m—1 \l—l / \^—1 / 

L 

= y]Tr(H,Hf) =LM (12) 

1=1 


Smaller values of / imply that the interfering subspaces 
are closely aligned and in the case of perfect alignment, 
/ = 0. In view of this, if the matrices correspond 

to the channels between a user and interfering base stations, 
/ can be defined as a measure of alignment for the user. 
The computation of the mean F and measure of alignment 
/ involve singular value decomposition (SVD) and hence 
are expensive to compute. In what follows, we will explore 
approximations for the measure of alignment function /. 





1) 3-Transmitter Case: In the case of 3-transmitter in¬ 
terference channels, each user has 2 interference subspaces 

2M 

and / has the form ^ Am(HiHf^ -f H 2 Hf^). In the 

m— M+1 

following Lemma, we will find closed form expression for 
the eigenvalues of HiHf^ -f H 2 H|^. 

Lemma 1. If Hi, H 2 G c 2 Mxm jj^g generator matrix of 
the subspaces Hi, 712 G G 2 m,m(C), eigenvalues of HiHf^ -f 
H2H|^ can be represented in descending order as 

1 -I- cos(0i), . . . , 1 -I- COs(0m), 1 — COs(0m), ..., 1 — cos(0i) 


that H'l and H 2 become ||3 


M 


M 


(13) 



( 14 ) 


(15) 


where 9m is the mth smallest principal angle between Hi and 

H2. 

L _ _ 

Proof: We will first show that eigenvalues of ^ 

Z —1 

are invariant under a common rotation transformation on the 
corresponding subspaces {Hi}l^i G Giv,m(C), i-e., H —>■ 
OH where O G is an arbitrary member of the 

orthogonal group 0{2M) IS). 

Let Hi ftp ^ V = Av and consider 


|^^(Oft,)(Oft,)^j v' = AV 
[^ftjftf ] O^v' = AV 


On right multiplying on both sides and using the fact that 
0^0 = Ijv, we have 


Thus ft'ifti^ -f ft 2 ft 2 ^ has the structure as in (fThl l. The 
eigenvalues of the matrix in (fThl) can be trivially found to be 
: 1 -I- cos(0i), 1 -I- cos(02), • ■ 1 + cos{9m), 1 — cos(0m), 

1 — cos(0M-i), • ■ 1 — cos(0i). This completes the proof. ■ 

It follows from Lemma [T| that 


2M 


M 


/= A„(HiHf+H 2 HP) = l-cos(0™) 

m=M+l m—1 

Since cos(^m) > cos^(^m) € [0;7r/2], we have that 

M 

^ 1 - cos(0^) < ^ 1 - cos^(6>m) 


(0^0) I^^H.ftf j O^v = A'O^v' 
> (^ft,ftf^ (O^v) = A'(O^V) 


Comparing with HiH^J v = Av, we have that A = 

A' and v = O^v'. 

Thus, it follows that 

A™ (ftiftf -f ft2ftP) = A„((Ofti)(Ofti)^ 

-G(Oft2)(Oft2)^) 

Let columns of fti and ft 2 be the corresponding principal 
vectors, and {hi}ff^, respectively. Define ft'i = 

Ofti and ft 2 = Oftp We can choose O G 0{2M) such 


M 

E 

m—1 m—1 

^/<d^(ftl,ft2) 

We can redefine / such that / = cig(fti,ft 2 ), which 
is intuitive as d^(ftl,ft2) is proportional to the degree of 
orthogonality between the corresponding subspaces and thus 
has the property that its value decreases as the interfering 
subspaces get closer or more aligned. 

2 ) General Case: Unlike the 3-transmitter case, it is diffi¬ 
cult to obtain closed-form expressions for eigenvalues in the 
general case where the network has K transmitters. In a iC 
transmitter network, there are {K—\) interference subspaces at 

2M /K-1 _ _ \ 

each user and thus / has the form ^ Am ( X) HjHp 

m—M+1 \ 1—1 

In the following Lemma, we extend the bound for the general 
case of iT-transmitter interference channels. 

Lemma 2. For iT-transmitter interference channels, 
the measure of alignment / is bounded above by 

. .min E dc(Hj,ft;) 

1^1 

Proof: 

Let F be the mean of Since F is the minimizer 

K-l 

of follows for any arbitrary I < j < K —1 

1^1 

that 


K-l 


K-l 


^ dliF,Ui) <Y d: 


1^1 


^^ft„fti 


/< 


K-l 

1^1 


(17) 













'1 + COS^(0i) 

0 

sin(26li)/2 

0 


Eh;h;^ = ■ 

0 

0 1-|-cos^(0m) 

0 

sin(20M)/2 

(16) 

sin(26»i)/2 

0 

sin^(0i) 

0 


0 

• • • sin(20M)/2) 

0 

sin^ (6 <m) 



Since ( [TtI i holds for any arbitrary j, it must hold for all 
j £ {1,..., itT — 1} and thus 

K-l 

(18) 

/-I 

This completes the proof. ■ 

K-l 

From (HB, we have f = dg(F,Hi), where F cor- 

1=1 

responds to the mean of the subspaces corresponding to 
{fii}f=i^- Since computing the mean F or even / (directly) 
is computationally prohibitive, we can approximate the mean 
F by an element in which is nearest to it. Indeed 

K-l 

H- where j = min argmin ^ ci^(Hj,H/) is closest to 

the mean F and we can redefine the measure of alignment as 
follows 

K-l 

/ = (19) 


follows 


K 


Tr E I U 

= Vm+1 (B„*) , Vm+2 (B„*) , . . . , V2M 


( 22 ) 


K 


where B^. = ^ Hjij. ;(H*. ;)^. The achievable sum- 

1=1,i^k 

rate 0 for the network is given by 


K 


Rsum — ^ ^ 


K 






k^l 


K 




1=1, l^k 


(23) 


V. Opportunistic Interference Alignment 
THROUGH User Pairing 


Note that this approximation to the actual measure of 
alignment is cheaper to compute. Also, for K = 3, the above 
expression reduces to the one obtained for the 3-transmitter 
interference channel. 


IV. Opportunistic Interference Alignment 
THROUGH User Selection 


In this section, we consider the user selection problem (refer 
Section Ill-Al l in which one user is selected in each cell such 
that opportunistic lA is achieved. The nth user in cell k 
calculates the measure of alignment function as follows 


K 

= min y] 

l<j<K 

j^k l—l,l^k 




( 20 ) 


where ; is an arbitrary generator matrix for ;. Fol¬ 

lowing this, each user feeds the measure of alignment back to 
its corresponding transmitter. After receiving this information 
from their users, transmitter k selects the user, n^, with the 
minimum value of measure of alignment 


nfe= argmin/* (21) 

l<n<K 


In this section, we consider the problem of finding 
transmitter-receiver pairing configuration (refer Section III-Bt 
in order to achieve lA opportunistically. For OIA in the user 
pairing framework, the receivers feed back the measure of 
alignment to a central node, which in turn decides the pairing 
configuration. 

Each user receives K, M-dimensional signals among which 
atmost one can be the desired signal. Unlike the OIA with 
user selection case, the desired signal is not predefined and it 
will depend on the channel conditions for all the users in the 
network. The measure of alignment at user n when it is paired 
with the fcth BS can be defined as 

K 

fn,k= min dl{Gn,j,Gn,i) (24) 

j^k l—l,l^k 

where Gn,i is an arbitrary generator matrix for Gn,i- 
Each user thus computes K measure of alignment functions, 
{fn,k}k=i corresponding to each BS. 

Let us define the vector of measure of alignment at user n, 
in as 

in = [fn,l, fn,2, . ■ . , fn,Kf (25) 


The selected user in cell k employs the post-processing 
matrix U*. which minimizes the interference leakage IB as 


Each user feedbacks its corresponding measure of alignment 
vectors {in}n=i to a central node. The center aggregates the 




















Fig. 3. 3-transmitter MIMO IC with M = 3 and N = 30 


data from all the users and forms the N x K feedback matrix, 
F defined as 

F = [fi, f2, f^]^ (26) 

Each entry in the matrix F corresponds to a pair in the 
original network. The smaller the value of the entry, the more 
likely it is for the corresponding link to have the interferences 
aligned and thus more likely to be chosen in the final user 
pairing solution. Having obtained the matrix F, the center can 
choose K non-conflicting pairs which constitute the minimum 
sum for the measure of alignment. Therefore, the optimization 
problem can be formulated as 


N K 


min 

p 

i=i j=i 

(27a) 

subject to 

j 

(27b) 


Y.p^,J = lyJ 

(27c) 



(27d) 


This optimization can be solved efficiently by the rectangular 
Hungarian algorithm IfTOl . After the optimal pairing configu¬ 
ration P* has been found, each user which is connected to a 
BS can employ a post-processing matrix which minimizes the 
interference leakage similar to (l22l l. Let be the user paired 
with BS k, i.e., P*. ^ = 1. The expression for achievable 
sum-rate will be same as (EDI. 

VI. Performance Comparison 

In this section, we compare the performance in terms 
of sum-rate and computational complexity of the proposed 
OIA algorithm with the conventional MAX-SNR and MIN- 
INR schemes 0, Q. MAX-SNR and MIN-INR have been 
proposed for the user selection framework in 0, Q. We 
extend them for the user pairing framework as done for OIA in 



Fig. 4. Complexity for 3-transmitter MIMO IC with M = 3 


Section In what follows, US and UP denote user selection 
and user pairing, respectively. 

A. Complexity Analysis 

In this section, we will discuss the computational complex¬ 
ity of the algorithms using flop counts. The complexity of an 
operation is counted as total number of flops required which is 
defined as a real floating point operation and we denote it by 
i/). The flop counts for some typical operations for a complex 


matrix A S ^j-g 

i/'(A -t- A) = 2mn (28a) 

iI){\\A\\f) = Amn (28b) 

V'(GSO(A)) = - 2mn (28c) 

i/;(SVD(A)) = 24mn2 -p -p (28d) 

'0(MUL(A)) = ip{AA^) = 8mn^ — 2mn (28e) 


where GSO stands for Gram-Schmidt orthogonalization, SVD 
stands for singular value decomposition and MUL(A) denotes 
the operation AA^. 

The proposed OIA in the user selection framework requires 
{K — 1) GSO operations, {K — 1) MUL operations and 
{^ 2 ^) matrix subtractions as well as ||.||ii’ operations at each 
user in every cell. At the selected user, MUL operation for 
{K — 1) times, {K — 2) matrix additions and a single SVD 
is required. Thus the total complexity for OIA in the user 
selection framework is given by 

iAoia-us = k(s {N%{4K - 4) + - llA + 8)) 

+ (A|(124 + 2A) + A|(A - 3)) ) (29) 

OIA in the user pairing framework requires K GSO op¬ 
erations, K MUL operations and (^) matrix subtractions as 
well as ||.||_F operations at each user. At the K selected users, 
MUL operation for {K — 1) times, {K — 2) matrix additions 






















and a single SVD is required. Thus the total complexity for 
OIA in the user pairing framework is given by 

V'OIA-UP = Nx (NliiAK) + - bK)) + 

K X (A^|(124 + 2K) + N]i{K - 3)) (30) 

The complexity for MAX-SNR and MIN-INR in the user 
selection framework denoted by V'max-snr-us and V'min-inr-us, 
respectively, is given in Q. MIN-INR in the user pairing 
framework requires K MUL operations, 2K matrix additions 
and a single SVD at every user. The total complexity for MIN- 
INR is given by 

V’MIN-INR-UP = Nx {Nl{\2SK) + NUiK)) (31) 

The MAX-SNR scheme requires K GSO operations and K 
SVD at every user in the user pairing framework. Thus, the 
total complexity for MAX-SNR is given by 

V'MAX-SNR-UP = Nx {Nl{12SK) - Nl{K)) (32) 

Note that we have ignored the complexity of solving the 
optimization problem in (l27T i which arises in the user pairing 
framework. This is because the computation happens only 
once at the center and not at the mobile users. 

Fig. 4 and Fig. 6 show the plot of computational complexity 
vs. total number of users N for a 3-transmitter MIMO IC 
with M — 3 and a 4-transmitter MIMO IC with M = 6, 
respectively. It can be observed that the complexity of OIA 
is only a small fraction of the complexity of MIN-INR and 
MAX-SNR schemes. Moreover, user pairing when compared 
to user selection has roughly the same complexity in case of 
proposed OIA, but the same is not true for both, MIN-INR 
and MAX-SNR. 

B. Sum-Rate 

Fig. 3 and Fig. 5 show the sum-rate vs. signal-to-noise 
ratio (SNR) plot for the proposed OIA and the conventional 
schemes for a 3-transmitter MIMO IC with M = 3 and 
a 4-transmitter MIMO IC with M = 6, respectively. It 
can be observed that the proposed OIA achieves sum-rates 
close to MIN-INR but at a significantly lower computational 
complexity. Moreover, with the same number of total users in 
the network, user pairing outperforms user selection. 

Fig. 7 shows the plot of sum-rate vs. the total number 
of users for a 4-transmitter MIMO IC with M = 6 and 
SNRs of 10 dB and 25 dB. As expected, the performance 
of all the algorithms improve as the number of users is 
increased. Also, user pairing provides more than a 4-fold 
gain over user selection in terms of total number of users 
required to achieve similar sum-rate performance. With respect 
to the number of users, the gaps in sum-rate performance for 
different algorithms is almost constant. Fig. 8 shows the plot 
of sum-rate vs. the number of antennas M for a 4-transmitter 
MIMO IC with N = 100 and SNRs of 10 dB and 25 dB. 
It can be observed that the sum-rate increases almost linearly 
with the number of transmit antennas for all the algorithms. 



Fig. 5. 4-transmitter MIMO IC with M = 6 and = 40 



Fig. 6. Complexity for 4-transmitter MIMO IC with M = 6 

VII. Conclusion 

In this paper, we have considered two different system mod¬ 
els, namely, user selection and user pairing for iT-transmitter 
MIMO interference channels. By exploiting multiuser diver¬ 
sity, we propose low complexity opportunistic interference 
alignment (OIA) algorithms for both the models. The proposed 
OIA algorithms are compared with conventional schemes, 
MIN-INR and MAX-SNR, and found to achieve comparable 
sum-rates but at a significantly reduced computational com¬ 
plexity. 
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